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Abstract 

Current theoretical approaches to the analysis of radiative heat exchange at the nanoscale are 
based on Rytov's stochastic electrodynamics. However, this approach falls short in the descrip- 
tion of microscale energy transfer since it overlooks non-resonant contributions arising from the 
coupling between different modes of relaxation in the material. We show that the phonon density 
of states given through a log-normal distribution accounts for such mode-coupling and leads to a 
general expression for the heat transfer coefficient which includes non-resonant contributions. This 
expression fits the existing experimental results with remarkable accuracy. Thus, our theory goes 
beyond stochastic electrodynamics and offers an overall explanation of energy transfer through 
micrometric gaps regardless of geometrical configurations. 
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1. Introduction. — Heat transfer between two objects at nanometric scales maintained at 
different temperatures was first studied following a stochastic or fluctuational electrodynam- 
ics formalism established by Rytov et al. pQ and Polder and van Hove [2 J . More recent studies 
have taken into account Casimir and van der Waals forces [3] , evanescent waves [U [5] , surface 
phonon polaritons j6] or have used Landauer-like reformulation [7j, to cite some examples. 
All these approaches have something in common: electrostatic interactions, linear response 
regime and, consequently, the fluctuation-dissipation theorem (FDT). Thus, an approach 
based on the contributions of fluctuating dipole effects seems to be the heart of a consid- 
erably simplified treatment of energy transfer at the nanoscale [21 El IS]- These approaches 
include a wide range of phenomena in which the energy transfer between molecules is dom- 
inated by dipole-dipole interactions, also known as Forster energy transfer [9j. Nonetheless, 
as two nanostructures, thermalized at different temperatures, come closer to each other, the 
distribution of charges and currents becomes asymmetric and therefore, defies description in 
terms of dipolar interactions. Hence, it is clear that one must bear in mind higher orders ef- 
fects beyond the dipole-dipole [ID] • In this sense, stochastic electrodynamics falls short in the 
whole description of microscale energy transfer since it overlooks non-resonant contributions 
to the system response which violates the FDT mentioned above [IT] . These non-resonant 
contributions come from the phonon-photon interaction and arise in the coupling between 
the different modes of relaxation in the material. In addition, some of the previous studies 
have apparently called into question the validity of the Derjaguin's approximation [51 IS] [T2] 
since an imbalance of action and reaction can occur in out-of-equilibrium situations [T5] . 

In this letter we present a mesoscopic treatment of this process in order to shed light 
on previous controversies. The assumption of a log-normal as the phonon density of states, 
which takes into account such mode-coupling, allows us to obtain a general expression for 
the heat-transfer coefficient including resonant and non-resonant contributions. Our theory 
goes beyond stochastic electrodynamics and offers an overall explanation of energy trans- 
fer through micrometric gaps regardless of geometrical configurations, fitting the existing 
experimental results with a high degree of accuracy. 

2. Mesoscopic non- equilibrium thermodynamics. — A thermodynamic description entails 
the formulation of the second law. This can be carried out by means of the Gibbs entropy 
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postulate [Hj , 

5(t) = ~ y n(p,t)ln^|dp + V (1) 

Here n(p,t) is the momentum distribution of photons, dp is the volume element in momen- 
tum space and V the volume. This gives us the non-equilibrium entropy of the photon gas 
plus the bath, where S eq is the equilibrium entropy and n eq the equilibrium distribution. 
Entropy is produced due to irreversible processes and its production is always positive, in 
accordance with the second law. Irreversible processes in non-equilibrium systems are de- 
scribed by means of currents, thermodynamic forces and the entropy production rate. The 
presence of a momentum distribution gradient dn(p,t)/dp yields a current J (p,t) which 
satisfies the relation J (p,t) = —Ddn(p,t)/dp, derived from the entropy production [Ti] . 
This relation accounts for the diffusion of probability between different regions of the phase 
space. Here, D(~ p 2 /time) is a kinetic coefficient, the diffusion coefficient in momentum 
space. 
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FIG. 1: (color online). Schematic model. Schematic diagram of the radiation exchanged between 
two materials maintained at different temperatures, T\ and T2, separated by a distance d. a, Two 
parallel surfaces, b, Sphere-surface system. 

This scheme can be applied to the description of the radiative heat transfer between two 



3 



surfaces at temperatures Ti (hot) and T 2 (cold), respectively, separated by a distance d, Fig. 
[TJi. Since in the gap between both surfaces there are only hot and cold photons interacting, 
it is plausible to assume that the average current in this gap results from the superposition 
of hot and cold currents, J Hit) and Jc(t), respectively, 

J(P, t) = Jc(t)5 (p - Pc ) + J H {t)5 (p - p H ) ■ (2) 

These currents are related to the net rate of emission of photons on both surfaces and thus 
must be proportional to the density of photons. When the characteristic length scales are 
comparable to the wavelength of thermal radiation, At = ch/kBT, the diffusion coefficient 
might depend on frequency. Here, h is the Planck's constant. For length scales d ^> Xt 
(i.e. low frequencies) the diffusion coefficient must be constant (the wave character disap- 
pears). On the other hand, when d < A^ (i.e. high frequencies) system-size effects become 
important and the diffusion coefficient should depend on the ratio Xx/d or equivalently on 
the frequency [TT], D(u) = (h 2 c/V) D(u). If the temperatures are kept constant, hot and 
cold photons will reach equilibrium with their respective baths. Hence, it follows from the 
integration of the previous general expression of the diffusion current with the help of Eq. 
(J2j) that a stationary value Jst(u) = D(uS) x [n(u, 2\) — n(ui,T 2 )], is established. Here, 
n(u,T)/V = 2N(u,T)/h 3 , with N(co,T) being the averaged number of quasi-particles in 
a elementary cell of volume h 3 of the phase-space given by the Planck's distribution [T5] , 
N(u,T) = 1/ [exp (fko/kT) - 1]. 

Since each photon carries an energy Huj, the heat flow Q follows from the sum of all the 
contributions as Q = (l/h 2 ) J HujJ s t(u})dp, where p = (fko/c) fl p . For the proposed case we 
assume that the density of vibrational states is achieved through the use of a log-normal 
distribution, which corresponds to 



JJ(u)) = -1= exp 

(Xu;ov27r 



ln(o;/u;o) 
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5(u-u R ), (3) 
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with k c = 1/Xt, a law which has been successful in a somewhat similar problem analyzed 
previously [TB]. It may be considered that the log-normal results from the existence of a 
hierarchy of relaxation mechanisms in the material [TT]. This law differs from the Debye 
approximation for the density of states u 2 /n 2 c 3 related to purely vibrational modes and 
is a characteristic of disordered systems to which dynamics is mainly due to slow relaxing 
modes, as for example glasses [TT] . 



Therefore, this contribution comes clearly from the excited phonons in the material due 
to the photon-phonon interaction. Here, u (characteristic frequency), and a (standard 
deviation) are two fitting parameters, and lor oc 1/d is a resonance frequency. In the 
blackbody radiation limit (far-field), D(u) = 1/4 for which we obtain Q = ob (Tf — ) 
which constitutes the Stefan-Boltzmann non-equilibrium law, with ob = 7f 2 k B /60h 3 c 2 being 
the Stefan constant. Nevertheless, there must be a multi-scale transition regime between the 
low-frequency and the high-frequency regimes for which Q contains contributions ^ n c n /d n , 
where c n corresponds to different length scales taking into account the material emissivity [3] , 
evanescent wave contributions [U E], Casimir and electrostatic forces [21 El E], and so on. 
This multi-scale regime is accounted for by assuming that D(oS) = J2 n X n uj n 5{(jj — ojr), 
where X n oc c n . As we have shown in Ref. [10] . the modification of the distribution of 
charges or currents due to proximity between objects induces multipolar interactions of 
an order higher than the dipolar interaction. This hierarchy of multipoles introduces a 
concomitant hierarchy of length scales. Accordingly, the multiscale regime mentioned above 
which includes the different contributions and effects referred to in the literature, clearly 
corresponds to a multipolar field expansion. 



In experiments, one usually measures the heat transfer coefficient H (ojr, Tq) = 
Q/ (T\ — T2), with Tq = (Ti + T2) /2, for which we obtain 



where k = k 2 kB/[(2Ti) 5 ^ 2 aujQ]. Here, we have taken into account terms up to the second 
order, involving far-field, electrostatic and out-of-equilibrium contributions. 

3. Results and discussion. — In order to test Eq. Q, we compare it with the near-field 
thermal radiation measurements between two parallel plates made of three different polar 
dielectric material combinations (S2O2 — SiC>2, SiC>2 — Si, and SiC>2 — Au) [6], [18]. Figure 
[2] shows a strong enhancement of the heat conductance for short distances, exceeding to the 
blackbody limit. Additionally, the contribution containing the log- normal in Eq. (|4]) has 
been plotted in all the shaded areas into the Fig. [2} 

Let us analyze the thermal conductance between a microsphere and a plate. In a previous 
discussion [6], the authors pointed out that Derjaguin's approximation [19] does not hold 



H (uJr, Tq) = AlObT'i 



(A + XiuJr + A 2 Wr) + 
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FIG. 2: (color online). Radiative heat transfer coefficients between two parallel plates as a function 
of the gap obtained by Shen et al. [6] . The black solid line is the limit of thermal radiation predicted 
by the usual blackbody radiation law, where the heat flux is calculated from Stefan-Boltzmann law 
at an average temperature T = 300 K. The black dashed lines show the analytical result obtained 
from Eq. Q by adjusting all parameters to the experimental data. The filled curves correspond 
to the non-Debye behavior for each of the material configurations, glass-glass, glass-silicon and 
glass-gold. 

for near- field radiative heat transfer. In contrast, other authors [3] demonstrated that this 
approach is in agreement with the data in the range 2.5 /im to 30 nm. Equation Q which 
avoids the use of Derjaguin's approximation, provides a general way which sheds light on 
the theory of radiative heat transfer at the nanoscale. Assuming the scheme in Fig. [TJo, 
the distance d between both surfaces depends on the curvature measured through the local 
radius r, as described in Ref. [5], ujr = Un(d(r)). Here, we use d — d+b + R— \/R 2 — r 2 , 
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where b represents a roughness parameter [5]. Hence, a more accurate measurement of the 
thermal conductance is given by the surface integral of Eq. (4), H = J Q R H (cur,T ) 2T\rdr. 
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FIG. 3: (color online). Thermal conductance between sphere and surface as function of the gap. 
a, Experimental data from the near-field radiation measurements between a silicon dioxide (glass) 
sphere and different substrate materials (glass, silicon, and gold) obtained by Shen et al. [Bj. The 
shaded areas correspond to non-Debye regime represented by the log-normal contribution in Eq. 
Q: blue for glass-glass, green for glass-silicon, and red for glass-gold material combinations, b, 
The red and green dots represent experimental data obtained by Rousseau et al. [5j with a sphere 
of different radius. 

From this, we calculate numerically the surface average by adjusting the parameters 
mentioned above to the experimental data obtained in the Refs. [5j E] . In Fig. [3^, we plotted 
the near- field conductance fitting Eq. (|4]) to the data [6]. Shaded areas reveal the influence 
of glassiness on the thermal conductance, which is quite similar to the contribution of non- 
resonant evanescent waves to near-field radiation [HE]. As a further check of the theory, 



we compare measurements with a microsphere of different radius [5], Fig. [3Jd. In the range 
30 nm to 2.5 y«m, the roughness parameter provides a slight deviation from 1/d asymptotic 
behavior, in such a way that the transition between the far and near-field occurs because of 
the non-Debye contribution, in contrast with previous conclusions [5]. Therefore, glassiness 
is a crucial ingredient in the near-field radiative heat exchange. 

Additional information about the value of the parameters for fitting Eq. Q to the 
experimental data can be found in the supplementary material. 

4- Conclusion. — In summary, we have evaluated thermal conductance in a wide range 
of length scales, from the far-field to the near-field, giving a thermokinetic description of 
several experiments involving heat radiation through a very narrow gap. Although near- 
field radiative transfer is a highly complex phenomenon, we have been able to provide a 
unified and highly accurate explanation of heat exchange processes at the nanoscale as well 
as a description of the transition between the far-field and the near-field regimes. Since the 
experiments examined may involve a great variety of nanostructures, our theory possesses 
a wide scope of applications. The general methodology presented here may also be used in 
the study of other heat exchange processes such as those occurring in phonon systems and 
in the analysis of thermal contributions to Casimir forces. 
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